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^5 ^ We analyse an implementation of a quantum computer using bosonic atoms in an 

optical lattice. We show that, even though the number of atoms per site and the 
G , tunneling rate between neighbouring sites is unknown, one may perform a universal 

^ ' set of gates by means of adiabatic passage. 
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1. Introduction 

> . 
^^ ^ Some tasks in Quantum Information require the implementation of quantum gates 

t:j- '• with a very high fidehty (A. M. Steane 2002; E. Knill 1998; D. Aharonov 1999). This 

implies that all parameters describing the physical system on which the computer 
is implemented have to be controlled with a very high precision, something which 
it is very hard to achieve in practice. 
f~^ ■ For example one can imagine an implementation in which qubits are stored in 

atoms and are manipulated using Raman transitions. It may happen that the rela- 
tive phase of the lasers driving a Raman transition can be controlled very precisely, 
whereas the corresponding Rabi frequency fJ has a larger uncertainty Af2. If we 
r^ . denote by T the time required to execute a local gate (of the order of ri~^), then 

d ' a high gate fidelity requires TAfi ^ 1 (equivalently, Afi/il ^ 1), which may very 

3 ,. hard to achieve, at least to reach the above mentioned threshold. 

In this paper we analyse an implementation of quantum computing using atoms 
confined in optical lattices. These systems have interesting features for quantum 
computing. Namely, a large number of atoms can be trapped in the lattice at a 
very low temperature, which provides a large number of qubits. Also, neutral atoms 
C^ \ interact weakly with the environment, which leads to a relatively slow decoherence. 

However, the same setups pose also important experimental challenges, such as 
being able to load the lattice with one atom per site or being able to measure 
the interaction and tunneling constants of these systems with high accuracy. These 
obstacles, together with the uncertainties in the atom-laser interaction, must be 
overcome to implement current proposals for quantum computing with neutral 
atoms (D. Jaksch 1999; R. lonicioiu 2002). 

In this paper we analyse a way that solves the above mentioned problems and 
show how to achieve a very high gate fidelity even when most of the parameters de- 
scribing the atomic ensemble (number of atoms per lattice site, tunneling rate, Rabi 
frequencies, etc) cannot be adjusted to precise values, and even have uncertainties 
of the order of the parameters themselves. 
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Our method combines the technique of adiabatic passage with ideas of quantum 
control theory. The use of adiabatic passage to implement quantum gates is not a 
new idea, and indeed several methods based on Berry phases have been proposed 
recently (P. Zanardi 1999; J. Pachos 2000; J. A. Jones 1999; G. Falci 2000; L.-M. 
Duan 2001). Furthermore, adiabatic passage techniques have been suggested as a 
way of implementing a universal set of holonomies (L.-M. Duan 2001), i.e. quantum 
gates which are carried out by varying certain parameters and whose outcome only 
depends on geometrical properties of the paths in parameter space (P. Zanardi 
1999; J. Pachos 2000). However, all these proposals are based on the existence of 
holonomies in the system, which in turn implies a huge degeneracy in the system. 
This will not be the case in our scheme. 

The outline of the paper is as follows. In §2 we introduce the requirements for 
quantum computing and show which tools are available in current experiments with 
cold atoms in optical lattices. We will demonstrate that, due to imperfections in 
the loading of the lattice, the Hamiltonian of the system is not known with enough 
accuracy to perform quantum computing in a 'traditional' way. In §3 we develop a 
technique to circumvent our ignorance about the Hamiltonian. Performing adiabatic 
passage with the different parameters of our problem, we show how to produce a 
universal set of gates (Hadamard, phase, and CNOT). In §4 we quantify the errors 
of our proposal, studying the influence of the speed of the adiabatic process, and 
of other imperfections. In §5 we summarize our results and offer some conclusions. 



2. Cold bosonic atoms in optical lattices 

We will consider a set of bosonic atoms confined in a periodic optical lattice at 
sufficiently low temperature such that only the first Bloch band is occupied. The 
atoms have two relevant internal (ground) levels, \a) and \h), and we wish to use 
this degree of freedom to store the qubit. This set-up has been studied in D. Jaksch 
(1999) where it has been shown how single quantum gates can be realized using 
lasers and two-qubit gates by displacing the atoms that are in a particular internal 
state to the next neighbour location. The basic ingredients of such a proposal have 
been recently realized experimentally (I. Bloch 2002, personal communication). 
However, in this and all other schemes so far (E. Charron 2002; K. Eckert 2002) 
it is assumed that there is a single atom per lattice site since otherwise even the 
concept of qubit is no longer valid. In present experiments, in which the optical 
lattice is loaded with a Bose-Einstein condensate (D. Jaksch 1998; M. Greincr 
2002), this only approximately true, since zero temperature is required and the 
number of atoms must be identical to the number of lattice sites. 



(a) Requirements for computation 

The uncertainty of the number of atoms per lattice site poses severe problems. 
Having ni atoms on the i-th cell, the configuration of this lattice size will be given 
by a combination of n^ + 1 possible states 

n^ = lini ^= ^ _ hfar^^vac),, k = Q...nA. (2.1) 
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To do quantum computing with m qubits, we must find a 2™-dimensional subspace 
He C Til (^ ■ ■ ■ <S) Tim, which is energetically separated from the rest, so that once 
we set our computer in a superposition of qubits |0) and |1), it docs not leave this 
subspace. A second, and stronger requirement is that our computation space must 
be an eigenspace of our Hamiltonian, with the same eigenvalue 

H{\zi)(g>---(g>\z„))^0, Vz, eZ2, (2.2) 

which we assume 0. Otherwise the trivial evolution of our system would spoil the 
quantum computation by introducing uncontrollable, unknown phases. 

Right from the beginning we forsee several difficulties. First, for arbitrary inter- 
actions, the states with different occupation numbers (2.1) will be regularly spaced 
and we will not be able to select our qubits. Furthermore, even if we customize the 
interactions between bosons, given that we basically ignore the number of atoms 
per site, n^, a basic requisite of our scheme will be to show that our procedure 
works independently of the occupation numbers. Both problems cannot be solved 
in general. We will rather have to impose some restrictions on our physical system, 
and this is the purpose of the following subsections. 

(b) Definition of qubit 

A crucial assumption which is suggested by the requirement (2.2), is to impose 
that the atoms in internal state \a) do not interact and do not hop to neighbouring 
sitesf . In the absence of external fields, the Hamiltonian describing our system is 



^-E 



-JbMi^i+ibk + blbk+i) + ■::Ubbb\b\bkbk 



(2.3) 



Here Ubh and J^jt describe the interactions between and the tunneling of atoms in 
state \b). We will assume that Jb,k can be set to zero and increased by adjusting 
the intensities of the lasers which create the optical lattice. 

For us a qubit will be formed by an aggregate of at least one atom per lattice 
site and the qubit basis will be formed by the states with at most one atom excited 
to the state b. More precisely, our computation will be performed in the space 



Tic{n^)--\\zl...ZM)^\{-^b''^^a^ ">ac) : Vz, G Z2 I . (2.4) 



i=l 



It is easy to check that for Jf, = 0, all our qubit states |5i . . . 6m) form degenerate 
linear eigenspace of our Hamiltonian (2.3), which is separated by an energy gap of 
Uhb from any other configuration. 

Definition (2.4) fulfills some of the requirements for quantum computing. How- 
ever, Tic does depend on the occupation number of the lattice, while a general state 
will be an incoherent superposition of different occupation numbers. It remains to 
show that we are able to produce quantum gates which are insensitive to the num- 
bers Ui. More precisely, if we design a protocol to produce the gate Uid, and this 
protocol is implemented by the unitary operation Ureal , we must prove that within 
the required accuracy 

t This may be achieved by tuning the scattering lengths and the optical lattice 
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The phases (fiijii) are irrelevant, since they are common to each of the possible 
computation spaces and final measurements will project our state to one of the 
subspaces Tic- 



(c) Available tools 

The quantum gates will be realized using lasers, switching the tunneling between 
neighbouring sites, and using the atom-atom interaction. We will now show how 
these elements introduce enough modifications to the Hamiltonian (2.3) so as to 
perform general quantum computing. 

For a single qubit gate on qubit k we can induce unitary transformations by 
means of Stark shifts and transitions between internal states. During the whole 
operation we set Jb.k = in order to isolate the qubits. The atom-laser interaction 
is then described by the Hamiltonian 

Hias,k = ^(4«fe - blbk) + ^(e"^46fe + e-'^lak). (2.6) 

For Ubb S> |Afe|, |rife|, we can replace (2.6) by an effective Hamiltonian of the form 

H, = ^a, + ^{a+e'^ + a_e-''^), (2.7) 

with A = Ak and f2 = flky/n^ — 1. 

For the realization of two-qubit operations we tilt the lattice using an electric 
field, Htiit = ^k^di'^k'^k + bkbk)- The tilting must be weak as to only virtual 
hopping of atoms of type \b) (Jb <C \Ubb — g\)- After adiabatic elimination we find 
that the effective Hamiltonian becomes 

i72 = A|ll)(ll|, A =—4^. (2.8) 

g- Ubb 

The Hamiltonians Hi and H2 pose now two problems. The first one is that Hi 
depends on the occupation numbers. A traditional approach to quantum computing 
would be to tune the parameters A and fi and let the resulting Hamiltonian operate 
for a time T, Ureal = exp{—iHiT). However, since the parameters are unknown, 
we cannot take this naive route. The second difficulty resides on the magnitude of 
Jb, g, and Ubb- these values are very sensitive to the properties of the lattice and 
difficult to control. At most, we will be able to assure that Jb and g are zero, or 
that g is similar to Ubb] but we will be unable to fix the value of A with enough 
accuracy that Ureal — exp{—iH2T) resembles a controUed-Z gate. 

In the following section we will solve these problems. In §4 develop an abstract 
protocol which, up from the Hamiltonians (2.7) and (2.8), produces a universal set 
of gates that can be used for quantum computing. Next in §4 we will study the 
influence of all the processes which we have neglected in the abstract derivation, 
such as interaction and hopping of atoms in state |a), sensitivity to occupation 
numbers, etc. 
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Figure 1. Schema of how the parameters of Hamihonian (2.7) have to be changed in 
order to perform a phase gate (a), and Hadamard gate (b). 



3. Computation with unknown parameters 

(a) Basic ideas 

Let us consider a set of qubits that can be manipulated according to the single 
qubit Hamiltonian (2.7) and the two-qubit Hamiltonian (2.8). We will assume that 
most of the parameters appearing in these Hamiltonians are basically unknown. On 
the other hand, wc will not consider any randomness in these parameters because 
the corresponding errors may be corrected with standard error correction methods 
(M. A. Nielsen 2002), as long as they are small, and in most cases random quick 
fluctuations of the parameters / will be averaged out in the adiabatic process. 

In particular wc will assume that only the phase of the laser, ip, can be precisely 
controlled. For the other parameters we will impose that: (i) they arc given by an 
unknown (single valued) function of some experimentally controllable parameters, 
(ii) they can be set to zero, and (iii) they are positivef . For example, we may have 
51 = f{I), where / is a parameter that can be experimentally controlled, and we 
only know about / that f{0) — and that we can reach some value Qm = .f{Im) 7^ 
for some Imt- Outside this, /(/) may change in different experimental realizations. 

The physical scenario described in §2 corresponds to this situation, but we want 
to stress that these conditions can be naturally met in more general scenarios. For 
example, the qubit states |0) and |1) may correspond to two degenerate atomic 
(ground state) levels which are driven by two lasers of the same frequency and 
different polarization. The corresponding Hamiltonian is given by (2.7), where the 
parameters if, fl, A describe the relative phase of the lasers, the Rabi frequency and 
detuning of the two-photon Raman transition, respectively. The Rabi frequency can 
be changed by adjusting the intensity of the lasers, and the detuning and the phase 
difference by using appropriate modulators. In practice, Q (A) can be set to zero 
very precisely by switching off the lasers (modulators) and (p may be very precisely 
controlled to any number between and 27r. However, fixing 17 or A to a precise 
value can be much more difficult. 

f This is just accommodate the physical restrictions of §2. The scheme actually becomes simpler 
when A or A may take negative values. 

I Note that, in many realistic implementations it is not possible to measure the dependence 
of these parameters (function /) because measurements are destructive (lead to heating or atom 
losses). Therefore /(/) is different in different experimental realizations. 
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The idea of obtaining perfect gates with unknown parameters combines the 
techniques of adiabatic passage (P. Zanardi 1999; J. Pachos 2000) with ideas of 
quantum control (L. Viola 1998; L.-M. Duan 1998). Let us briefly recall the adi- 
abatic theorem, which is a fundamental tool in our method. Suppose we have a 
Hamiltonian that depends parametrically on a set of parameters, denoted by p, 
which are changed adiabatically with time along a given trajectory p{t). After a 
time T, the unitary operator corresponding to the evolution is 

[/(T) =^e'(*°+'^-)|<I>„b(r)])(<I>„[MO)]|. (3.1) 

Here, \^a{'p)) a-re the eigenstates of the Hamiltonian for which the parameters take 
on the values p. The phase 4>a is a dynamical phase that explicitly depends on how 
the parameters p are changed with time, whereas the phase ^„ is a purely geomet- 
rical phase ands depends on the trajectory described in the parameter space. Our 
basic idea to perform any given gate is first to design the change of the param- 
eters in the Hamiltonians (2.7)-(2.8) such that the eigenvectors evolve according 
to the desired gate, and then to repeat the procedure changing the parameters 
appropriately in order to cancel the geometric and dynamical phases. 

[b) Local gates 

Using the previous ideas we are able to implement a universal set of gates, 
which is made of a phase gate, U ~ e**'^^/^, a Hadamard gate and a CNOT gate. 
To perform the phase gate U ~ g^s<^i/2 -^q work with the single-qubit Hamiltonian 
(2.7). We set A = for all times and change the remaining parameters {^,(p) as 
depictured in figure 1(a): 

(0, 0) ^ (r!™, 0) ^^ (a„, e/2) ^^ (f7„, e/2 + n) 

^^^ {n^,9 + Tr)^^{0,9 + 7T) (3.2) 

All steps are performed adiabatically and require a total time T, except for step 
(iii) whose double arrow indicates a sudden change of parameters. Note that J1(0) = 
n{2T) = 0, n{t) = 0(2T - t) and (p{t) ^ tt + - ip{2T - i), which does not 
require the knowledge of the function / but implies a precise control of the phase. 
A simple analysis shows that (i-v) achieve the desired transformation |0) -^ e*^/^|0), 
|1) -^ e^*^/^|l). Note also that the dynamical and geometrical phases acquired in 
the adiabatic processes (i-v) cancel out. 

The Hadamard gate can be performed in a similar way. In the space [A, ilj; = 
ncos((/5)], the protocol is 

(0,a„) ^ (A„,O„)*-^^(A,„,0)*^^A„„l]„) 

<^' (0, a^) ^ (0, -n^) '^^ (A, -o„) <^^^ (A, 0), (3.3) 

as shown in figure l(b-c). In order to avoid the dynamical phases, we have to make 
sure that steps (i-v) are run in half the time as (vi-vii). More precisely, if t < T, 
we must ensure that A(t) = A(T - i), 17^(i) = n^{T - t), A(T + t) = A(i/2) and 
n^{T + t) = r2^(t/2). With this requisite we get -^{\0) + \l))^ |0), ■j^{\0)-\l)) ^ 
— 11). Again, the whole procedure does not require us to know Q or A, but rather 
to control the evolution of the experimental parameters which determine them. 
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Figure 2. Log-log plot of the gate error, E = 1 — JT, for the Hadamard 
(solid), phase (dashed) and CNOT (dotted) gates, when we use the ideal Hamil- 
tonians (2.7) and (2.8) to simulate them. The parameters for the simulations are 
{Am/ilm = Am/fim — 1/10, ifm = 7r/4}, aud we chauge the speed of the process to 
study the influence of non-adiabaticity. 

(c) Nonlocal gates 

The C-NOT gate requires the combination of two two-qubit processes using H2 
and one local gate. The first process involves changing the parameters [A, Q^ = 
ricos((^)] of equation (2.8) according to 



(Am,0) 



(iv) 



^^{o,n^)^U\o,-n^) 






0). 



(3.4) 



(A„,-f7 
This procedure gives rise to the transformation 

Ui^\0){0\(g)l+e'i\l){l\(g)icry, (3.5) 

where C = /o 5{t)dt is an unknown dynamical phase. The second operation required 
is a NOT on the^ first qubit U2 = (|0)(1| + |1)(0|) (g)I. Finally, if A(i)(i) denotes 
the evolution of A in equation (3.4), we need to follow a path such that A('^)(i) = 
A(i)(i), n(3)(^) ^ 0. If the timing is correct, we achieve C/g = (|0)(0| + e*«|l)(l|) ®1 
Everything combined gives us the CNOT up to a global unimportant phase Ucnot — 
|0)(0| 01+ |1)(1| O iay = e-'iU2U3U2Ui. 

4. Control of errors 

In this section we study whether it is feasible to apply the methods developed 
in §3 to the physical setup envisioned in §2. First of all we will study how fast 
the operations from §3 must be performed in order to minimize the deviations 
from the adiabatic theorem. And second and most important, we have to consider 
contributions to the energy which escape the terms considered in equations (2.7) 
and (2.8). We will analyse both sources of error separately, combining analytical 
estimates with numerical simulations of our techniques for small number of atoms. 

(a) Adiabaticity 

To study the sensitivity of our method against non-adiabatic processes, we have 
simulated the protocols (3.2)-(3.4) using the ideal model given by Hamiltonians 
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Figure 3. Log-log plot of the gate error, E = 1 — JT, for the (a) phase, (b) Hadamard 
and (c) CNOT gates. The parameters for the simulations are Uaa = Uab = Jm, 
Ja = Jb,g = Ubb + Uab/'i, and T — 100/f2m- For the local gates we choose Am = 60^ ~ 1, 
and for the nonlocal gate Qm ~ Jm/^- For each simulation we choose different population 
imbalance (|n — 7n| =0,i,2 for solid, dashed and dotted lines), and change the interaction 
constant Ubb- 

(2.7) and (2.8). For each of the gates we have fixed all parameters {Am/^m = 
Am/^m = 1/10, ipm = 7r/4} except the time, and then we have computed how 
the error decreases as we decrease the speed of the adiabatic passage. The results 
are shown in figure 2. As a figure of merit wc have chosen the gate fidelity (M. A. 
Nielsen 2002) 

(4.1) 



:F^2-\TT{Ul^^,Ureal}? 



where n is the number of qubits involved in the gate, Uideai is the gate that we wish 
to produce and Ureal is the actual operation performed. As expected, the adiabatic 
theorem applies when the processes are performed with a sufficiently slow speed. 
Typically a time T ^ 300/ilm is required for the desired fidelity !F = 1 — ID"**. 

It is also worth mentioning, that in figure 2 and 3, strong, rapid oscillations of the 
error are seen. These oscillations are due to cither the nonadiabaticity of the process 
(figure 2), or to imperfections in the Hamiltonian (figure 3). For a two-level system 
undergoing adiabatic evolution it is easy to prove that, while the amplitude of the 
oscillations is proportional to the speed of the adiabatic process, the frequency 
is instead related to the energy difference between contiguous eigenspaces. This 
frequency is for us unknown, and consequently, these oscillations may not be used 
to improve the accuracy of our method by looking for some 'magic times'. 

(&) Imperfections 

Outside the non-adiabaticity of a real experiment, there are two other sources of 
error which we must consider, (i) The quotient Jt/Ubb is nonzero, which means that 
more than one atom per well may be excited, (ii) Atoms in state \a) interact and hop. 
This introduces new terms in equation (2.3), which are of the form Ja, fc(ajj,afc+i + 
O'k+i'^k), UaaO-k'^k'^kO-k, and UabO^bkO-kbk- And finally, (iii) atoms in either state may 
jump to neighbouring sites, permanently changing the occupation numbers. 

The effects (i) and (iii) are suppressed if the coupling between internal levels 
and the amplitude of hopping are both small compared to the energy gap between 
our computation space and the undcsired excitations. In other words, we require 



{n/Ubbf<^l and [Ja^k/Ubb? <^l- 



(4.2) 
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To analyse the remaining errors we develop an effective Hamiltonian which contains 
(2.6) and (2.3) plus new terms {Ja,Uab,Uaa) that we did not consider before. In 
equation (2.6), the virtual excitation of two atoms increments the parameter A by 
an unknown amount, Aoff ~ A + 2ft'^nk/{A + Uab ~ Ubb)- If 

Uab < Ubb and n^UkT/Ubb < 1, (4.3) 

this shift may be neglected. In the two-qubit gates the energy shifts are instead 
due to virtual hopping of all types of atoms, and they are also accompanied by the 
possibility of swapping both qubits (|01) ^^ |10))- Both contributions arc of the 
order of max(J6, Ja)^ jg^ ^ ,P jUbb, and for 

j'^TjUbb « 1 (4.4) 

they also may be neglected. 

To quantitatively determine the influence of these errors we have simulated 
the evolution of two atomic ensembles with an effective Hamiltonian which results 
of applying second order perturbation theory to equation (2.3), and which takes 
into account all important processes. The results are shown in figures 3. In these 
pictures we show the error of the gates for simulations in which all parameters are 
fixed, except for Ubb and the occupation numbers of the wells. The first conclusion 
is that the stronger the interaction between atoms in state |6), the smaller the 
energy shifts. This was already evident from our analytical estimates, because all 
errors are proportional to 1/Ubb- Typically, a ratio Ubb — 10'*C/afc is required to 
make T = \ — 10^''. Second, the larger the number of atoms per well, the poorer 
the fidelity of the local gates [Figure 3(a-b)]. And finally, as figure 3(c) shows, the 
population imbalance between wells influences very little the fidelity of the two- 
qubit gate. 

5. Conclusions 

In this work we have shown that it is possible to perform quantum computation 
with cold atoms in a tunable optical lattice. Our scheme is based on performing 
adiabatic passage with one-qubit (2.7) and two-qubit (2.8) Hamiltonians. With 
selected paths and appropriate timing, it is possible to perform a universal set of 
gates (Hadamard gate, phase gate and a CNOT). Thanks to the adiabatic passage, 
the proposal works even when the number of atoms per lattice is unknown or the 
constants in the governing Hamiltonians have large uncertainties. These procedures 
can not only be used for quantum computing but also for quantum simulation (E. 
Jane 2002), and the same ideas can also be applied to other setups like the micro- 
traps demonstrated in R. Dumke (2002). 

We thank D. Liebfried and P. ZoUer for discussions and the EU project EQUIP (contract 
IST-1999-11053). 
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